SARS-CoV-2 spike protein diversity at an intra-host level, among SARS-CoV-2 infected individuals in South Africa, 2020 to 2022.

Intra-host diversity studies are used to characterise the mutational heterogeneity of SARS-CoV-2 infections in order to understand the impact of virus-host adaptations. This study investigated the frequency and diversity of the spike (S) protein mutations within SARS-CoV-2 infected South African individuals. The study included SARS-CoV-2 respiratory samples, from individuals of all ages, received at the National Health Laboratory Service at Charlotte Maxeke Johannesburg Academic hospital, Gauteng, South Africa, from June 2020 to May 2022. Single nucleotide polymorphism (SNP) assays and whole genome sequencing were performed on a random selection of SARS-CoV-2 positive samples. The allele frequency (AF) was determined using TaqMan Genotyper software for SNP PCR analysis and galaxy.eu for analysis of FASTQ reads from sequencing. The SNP assays identified 5.3% (50/948) of Delta cases with heterogeneity at delY144 (4%; 2/50), E484Q (6%; 3/50), N501Y (2%; 1/50) and P681H (88%; 44/50), however only heterogeneity for E484Q and delY144 were confirmed by sequencing. From sequencing we identified 9% (210/2381) of cases with Beta, Delta, Omicron BA.1, BA.2.15, and BA.4 lineages that had heterogeneity in the S protein. Heterogeneity was primarily identified at positions 19 (1.4%) with T19IR (AF 0.2-0.7), 371 (92.3%) with S371FP (AF 0.1-1.0), and 484 (1.9%) with E484AK (0.2-0.7), E484AQ (AF 0.4-0.5) and E484KQ (AF 0.1-0.4). Mutations at heterozygous amino acid positions 19, 371 and 484 are known antibody escape mutations, however the impact of the combination of multiple substitutions identified at the same position is unknown. Therefore, we hypothesise that intra-host SARS-CoV-2 quasispecies with heterogeneity in the S protein facilitate competitive advantage of variants that can completely/partially evade host's natural and vaccine-induced immune responses.

Molecular Epidemiology of SARS-CoV-2 during Five COVID-19 Waves and the Significance of Low-Frequency Lineages.

SARS-CoV-2 lineages and variants of concern (VOC) have gained more efficient transmission and immune evasion properties with time. We describe the circulation of VOCs in South Africa and the potential role of low-frequency lineages on the emergence of future lineages. Whole genome sequencing was performed on SARS-CoV-2 samples from South Africa. Sequences were analysed with Nextstrain pangolin tools and Stanford University Coronavirus Antiviral & Resistance Database. In 2020, 24 lineages were detected, with B.1 (3%; 8/278), B.1.1 (16%; 45/278), B.1.1.348 (3%; 8/278), B.1.1.52 (5%; 13/278), C.1 (13%; 37/278) and C.2 (2%; 6/278) circulating during the first wave. Beta emerged late in 2020, dominating the second wave of infection. B.1 and B.1.1 continued to circulate at low frequencies in 2021 and B.1.1 re-emerged in 2022. Beta was outcompeted by Delta in 2021, which was thereafter outcompeted by Omicron sub-lineages during the 4th and 5th waves in 2022. Several significant mutations identified in VOCs were also detected in low-frequency lineages, including S68F (E protein); I82T (M protein); P13L, R203K and G204R/K (N protein); R126S (ORF3a); P323L (RdRp); and N501Y, E484K, D614G, H655Y and N679K (S protein). Low-frequency variants, together with VOCs circulating, may lead to convergence and the emergence of future lineages that may increase transmissibility, infectivity and escape vaccine-induced or natural host immunity.

Performance of rapid antigen tests in identifying Omicron BA.4 and BA.5 infections in South Africa.

BACKGROUND: Concerns around accuracy and performance of rapid antigen tests continue to be raised with the emergence of new SARS-CoV-2 variants. OBJECTIVE: To evaluate the performance of two widely used SARS-CoV-2 rapid antigen tests during BA.4/BA.5 SARS-CoV-2 wave in South Africa (May - June 2022). STUDY DESIGN: A prospective field evaluation compared the SARS-CoV-2 Antigen Rapid test from Hangzhou AllTest Biotech (nasal swab) and the Standard Q COVID-19 Rapid Antigen test from SD Biosensor (nasopharyngeal swab) to the Abbott RealTime SARS-CoV-2 assay (nasopharyngeal swab) on samples collected from 540 study participants. RESULTS: Overall 28.52% (154/540) were SARS-CoV-2 RT-PCR positive with median cycle number value of 12.30 (IQR 9.30-19.40). Out of the 99 successfully sequenced SARS-CoV-2 positive samples, 18 were classified as BA.4 and 56 were classified as BA.5. The overall sensitivities of the AllTest SARS-CoV-2 Ag test and Standard Q COVID-19 Ag test were 73.38% (95% CI 65.89-79.73) and 74.03% (95% CI 66.58-80.31) and their specificities were 97.41% (95% CI 95.30-98.59) and 99.22% (95% CI 97.74-99.74) respectively. Sensitivity was >90% when the cycle number value was <20. The sensitivity of both rapid tests was >90% in samples infected with Omicron sub-lineage BA.4 and BA.5. CONCLUSION: Accuracy of tested rapid antigen tests that target the nucleocapsid SARS-CoV-2 protein, were not adversely affected by BA.4 and BA.5 Omicron sub-variants.

Evaluation of miniaturized Illumina DNA preparation protocols for SARS-CoV-2 whole genome sequencing.

The global pandemic caused by SARS-CoV-2 has increased the demand for scalable sequencing and diagnostic methods, especially for genomic surveillance. Although next-generation sequencing has enabled large-scale genomic surveillance, the ability to sequence SARS-CoV-2 in some settings has been limited by the cost of sequencing kits and the time-consuming preparations of sequencing libraries. We compared the sequencing outcomes, cost and turn-around times obtained using the standard Illumina DNA Prep kit protocol to three modified protocols with fewer clean-up steps and different reagent volumes (full volume, half volume, one-tenth volume). We processed a single run of 47 samples under each protocol and compared the yield and mean sequence coverage. The sequencing success rate and quality for the four different reactions were as follows: the full reaction was 98.2%, the one-tenth reaction was 98.0%, the full rapid reaction was 97.5% and the half-reaction, was 97.1%. As a result, uniformity of sequence quality indicated that libraries were not affected by the change in protocol. The cost of sequencing was reduced approximately seven-fold and the time taken to prepare the library was reduced from 6.5 hours to 3 hours. The sequencing results obtained using the miniaturised volumes showed comparability to the results obtained using full volumes as described by the manufacturer. The adaptation of the protocol represents a lower-cost, streamlined approach for SARS-CoV-2 sequencing, which can be used to produce genomic data quickly and more affordably, especially in resource-constrained settings.

Genomic epidemiology of SARS-CoV-2 during the first four waves in Mozambique

Mozambique reported the first case of coronavirus disease 2019 (COVID-19) in March 2020 and it has since spread to all provinces in the country. To investigate the introductions and spread of SARS-CoV-2 in Mozambique, 1 142 whole genome sequences sampled within Mozambique were phylogenetically analyzed against a globally representative set, reflecting the first 25 months of the epidemic. The epidemic in the country was marked by four waves of infection, the first associated with B.1 ancestral lineages, while the Beta, Delta, and Omicron Variants of Concern (VOCs) were responsible for most infections and deaths during the second, third, and fourth waves. Large-scale viral exchanges occurred during the latter three waves and were largely attributed to southern African origins. Not only did the country remain vulnerable to the introductions of new variants but these variants continued to evolve within the borders of the country. Due to the Mozambican health system already under constraint, and paucity of data in Mozambique, there is a need to continue to strengthen and support genomic surveillance in the country as VOCs and Variants of interests (VOIs) are often reported from the southern African region.

Molecular Epidemiology and Diversity of SARS-CoV-2 in Ethiopia, 2020-2022.

Ethiopia is the second most populous country in Africa and the sixth most affected by COVID-19 on the continent. Despite having experienced five infection waves, >499,000 cases, and ~7500 COVID-19-related deaths as of January 2023, there is still no detailed genomic epidemiological report on the introduction and spread of SARS-CoV-2 in Ethiopia. In this study, we reconstructed and elucidated the COVID-19 epidemic dynamics. Specifically, we investigated the introduction, local transmission, ongoing evolution, and spread of SARS-CoV-2 during the first four infection waves using 353 high-quality near-whole genomes sampled in Ethiopia. Our results show that whereas viral introductions seeded the first wave, subsequent waves were seeded by local transmission. The B.1.480 lineage emerged in the first wave and notably remained in circulation even after the emergence of the Alpha variant. The B.1.480 was outcompeted by the Delta variant. Notably, Ethiopia's lack of local sequencing capacity was further limited by sporadic, uneven, and insufficient sampling that limited the incorporation of genomic epidemiology in the epidemic public health response in Ethiopia. These results highlight Ethiopia's role in SARS-CoV-2 dissemination and the urgent need for balanced, near-real-time genomic sequencing.

Clinical Evaluation of Severe Acute Respiratory Syndrome Coronavirus 2 Rapid Antigen Tests During the Omicron Wave in South Africa.

We evaluated the performance of nasal and nasopharyngeal Standard Q COVID-19 [coronavirus disease 2019] Ag tests (SD Biosensor) and the Panbio COVID-19 Ag Rapid Test Device (nasal; Abbott) against the Abbott RealTime severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) assay during the Omicron (clades 21M, 21K, and 21L) wave in South Africa. Overall, all evaluated tests performed well, with high sensitivity (range, 77.78%-81.42%) and excellent specificity values (>99%). The sensitivity of rapid antigen tests increased above 90% in samples with cycle threshold <20, and all 3 tests performed best within the first week after symptom onset.

Omicron BA.4/BA.5 escape neutralizing immunity elicited by BA.1 infection.

SARS-CoV-2 Omicron (B.1.1.529) BA.4 and BA.5 sub-lineages, first detected in South Africa, have changes relative to Omicron BA.1 including substitutions in the spike receptor binding domain. Here we isolated live BA.4 and BA.5 viruses and measured BA.4/BA.5 neutralization elicited by BA.1 infection either in the absence or presence of previous vaccination as well as from vaccination without BA.1 infection. In BA.1-infected unvaccinated individuals, neutralization relative to BA.1 declines 7.6-fold for BA.4 and 7.5-fold for BA.5. In vaccinated individuals with subsequent BA.1 infection, neutralization relative to BA.1 decreases 3.2-fold for BA.4 and 2.6-fold for BA.5. The fold-drop versus ancestral virus neutralization in this group is 4.0-fold for BA.1, 12.9-fold for BA.4, and 10.3-fold for BA.5. In contrast, BA.4/BA.5 escape is similar to BA.1 in the absence of BA.1 elicited immunity: fold-drop relative to ancestral virus neutralization is 19.8-fold for BA.1, 19.6-fold for BA.4, and 20.9-fold for BA.5. These results show considerable escape of BA.4/BA.5 from BA.1 elicited immunity which is moderated with vaccination and may indicate that BA.4/BA.5 may have the strongest selective advantage in evading neutralization relative to BA.1 in unvaccinated, BA.1 infected individuals.

Emergence of SARS-CoV-2 Omicron lineages BA.4 and BA.5 in South Africa.

Three lineages (BA.1, BA.2 and BA.3) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant of concern predominantly drove South Africa's fourth Coronavirus Disease 2019 (COVID-19) wave. We have now identified two new lineages, BA.4 and BA.5, responsible for a fifth wave of infections. The spike proteins of BA.4 and BA.5 are identical, and similar to BA.2 except for the addition of 69-70 deletion (present in the Alpha variant and the BA.1 lineage), L452R (present in the Delta variant), F486V and the wild-type amino acid at Q493. The two lineages differ only outside of the spike region. The 69-70 deletion in spike allows these lineages to be identified by the proxy marker of S-gene target failure, on the background of variants not possessing this feature. BA.4 and BA.5 have rapidly replaced BA.2, reaching more than 50% of sequenced cases in South Africa by the first week of April 2022. Using a multinomial logistic regression model, we estimated growth advantages for BA.4 and BA.5 of 0.08 (95% confidence interval (CI): 0.08-0.09) and 0.10 (95% CI: 0.09-0.11) per day, respectively, over BA.2 in South Africa. The continued discovery of genetically diverse Omicron lineages points to the hypothesis that a discrete reservoir, such as human chronic infections and/or animal hosts, is potentially contributing to further evolution and dispersal of the virus.

Identification of SARS-CoV-2 Omicron variant using spike gene target failure and genotyping assays, Gauteng, South Africa, 2021.

The circulation of Omicron BA.1 led to the rapid increase in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cases in South Africa in November 2021, which warranted the use of more rapid detection methods. We, therefore, assessed the ability to detect Omicron BA.1 using genotyping assays to identify specific mutations in SARS-CoV-2 positive samples, Gauteng province, South Africa. The TaqPath™ COVID-19 real-time polymerase chain reaction assay was performed on all samples selected to identify spike gene target failure (SGTF). SARS-CoV-2 genotyping assays were used for the detection of del69/70 and K417N mutation. Whole-genome sequencing was performed on a subset of genotyped samples to confirm these findings. Of the positive samples received, 11.0% (175/1589) were randomly selected to assess if SGTF and genotyping assays, that detect del69/70 and K417N mutations, could identify Omicron BA.1. We identified SGTF in 98.9% (173/175) of samples, of which 88.0% (154/175) had both the del69/70 and K417N mutation. The genotyped samples (45.7%; 80/175) that were sequenced confirmed Omicron BA.1 (97.5%; 78/80). Our data show that genotyping for the detection of the del69/70 and K417N coupled with SGTF is efficient to exclude Alpha and Beta variants and rapidly detect Omicron BA.1. However, we still require assays for the detection of unique mutations that will allow for the differentiation between other Omicron sublineages. Therefore, the use of genotyping assays to detect new dominant or emerging lineages of SARS-CoV-2 will be beneficial in limited-resource settings.

Omicron infection enhances Delta antibody immunity in vaccinated persons.

The extent to which Omicron infection1-9, with or without previous vaccination, elicits protection against the previously dominant Delta (B.1.617.2) variant is unclear. Here we measured the neutralization capacity against variants of severe acute respiratory syndrome coronavirus 2 in 39 individuals in South Africa infected with the Omicron sublineage BA.1 starting at a median of 6 (interquartile range 3-9) days post symptom onset and continuing until last follow-up sample available, a median of 23 (interquartile range 19-27) days post symptoms to allow BA.1-elicited neutralizing immunity time to develop. Fifteen participants were vaccinated with Pfizer's BNT162b2 or Johnson & Johnson's Ad26.CoV2.S and had BA.1 breakthrough infections, and 24 were unvaccinated. BA.1 neutralization increased from a geometric mean 50% focus reduction neutralization test titre of 42 at enrolment to 575 at the last follow-up time point (13.6-fold) in vaccinated participants and from 46 to 272 (6.0-fold) in unvaccinated participants. Delta virus neutralization also increased, from 192 to 1,091 (5.7-fold) in vaccinated participants and from 28 to 91 (3.0-fold) in unvaccinated participants. At the last time point, unvaccinated individuals infected with BA.1 had low absolute levels of neutralization for the non-BA.1 viruses and 2.2-fold lower BA.1 neutralization, 12.0-fold lower Delta neutralization, 9.6-fold lower Beta variant neutralization, 17.9-fold lower ancestral virus neutralization and 4.8-fold lower Omicron sublineage BA.2 neutralization relative to vaccinated individuals infected with BA.1. These results indicate that hybrid immunity formed by vaccination and Omicron BA.1 infection should be protective against Delta and other variants. By contrast, infection with Omicron BA.1 alone offers limited cross-protection despite moderate enhancement.

Comparison of SARS-CoV-2 sequencing using the ONT GridION and the Illumina MiSeq.

BACKGROUND: Over 4 million SARS-CoV-2 genomes have been sequenced globally in the past 2 years. This has been crucial in elucidating transmission chains within communities, the development of new diagnostic methods, vaccines, and antivirals. Although several sequencing technologies have been employed, Illumina and Oxford Nanopore remain the two most commonly used platforms. The sequence quality between these two platforms warrants a comparison of the genomes produced by the two technologies. Here, we compared the SARS-CoV-2 consensus genomes obtained from the Oxford Nanopore Technology GridION and the Illumina MiSeq for 28 sequencing runs. RESULTS: Our results show that the MiSeq had a significantly higher number of consensus genomes classified by Nextclade as good and mediocre compared to the GridION. The MiSeq also had a significantly higher genome coverage and mutation counts than the GridION. CONCLUSION: Due to the low genome coverage, high number of indels, and sensitivity to SARS-CoV-2 viral load noted with the GridION when compared to MiSeq, we can conclude that the MiSeq is more favourable for SARS-CoV-2 genomic surveillance, as successful genomic surveillance is dependent on high quality, near-whole consensus genomes.

Emergence and phenotypic characterization of the global SARS-CoV-2 C.1.2 lineage.

Global genomic surveillance of SARS-CoV-2 has identified variants associated with increased transmissibility, neutralization resistance and disease severity. Here we report the emergence of the PANGO lineage C.1.2, detected at low prevalence in South Africa and eleven other countries. The initial C.1.2 detection is associated with a high substitution rate, and includes changes within the spike protein that have been associated with increased transmissibility or reduced neutralization sensitivity in SARS-CoV-2 variants of concern or variants of interest. Like Beta and Delta, C.1.2 shows significantly reduced neutralization sensitivity to plasma from vaccinees and individuals infected with the ancestral D614G virus. In contrast, convalescent donors infected with either Beta or Delta show high plasma neutralization against C.1.2. These functional data suggest that vaccine efficacy against C.1.2 will be equivalent to Beta and Delta, and that prior infection with either Beta or Delta will likely offer protection against C.1.2.

Rapid Replacement of SARS-CoV-2 Variants by Delta and Subsequent Arrival of Omicron, Uganda, 2021.

Genomic surveillance in Uganda showed rapid replacement of severe acute respiratory syndrome coronavirus 2 over time by variants, dominated by Delta. However, detection of the more transmissible Omicron variant among travelers and increasing community transmission highlight the need for near-real-time genomic surveillance and adherence to infection control measures to prevent future pandemic waves.

Selection Analysis Identifies Clusters of Unusual Mutational Changes in Omicron Lineage BA.1 That Likely Impact Spike Function.

Among the 30 nonsynonymous nucleotide substitutions in the Omicron S-gene are 13 that have only rarely been seen in other SARS-CoV-2 sequences. These mutations cluster within three functionally important regions of the S-gene at sites that will likely impact (1) interactions between subunits of the Spike trimer and the predisposition of subunits to shift from down to up configurations, (2) interactions of Spike with ACE2 receptors, and (3) the priming of Spike for membrane fusion. We show here that, based on both the rarity of these 13 mutations in intrapatient sequencing reads and patterns of selection at the codon sites where the mutations occur in SARS-CoV-2 and related sarbecoviruses, prior to the emergence of Omicron the mutations would have been predicted to decrease the fitness of any virus within which they occurred. We further propose that the mutations in each of the three clusters therefore cooperatively interact to both mitigate their individual fitness costs, and, in combination with other mutations, adaptively alter the function of Spike. Given the evident epidemic growth advantages of Omicron overall previously known SARS-CoV-2 lineages, it is crucial to determine both how such complex and highly adaptive mutation constellations were assembled within the Omicron S-gene, and why, despite unprecedented global genomic surveillance efforts, the early stages of this assembly process went completely undetected.

Targeted Sanger sequencing to recover key mutations in SARS-CoV-2 variant genome assemblies produced by next-generation sequencing.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is adaptively evolving to ensure its persistence within human hosts. It is therefore necessary to continuously monitor the emergence and prevalence of novel variants that arise. Importantly, some mutations have been associated with both molecular diagnostic failures and reduced or abrogated next-generation sequencing (NGS) read coverage in some genomic regions. Such impacts are particularly problematic when they occur in genomic regions such as those that encode the spike (S) protein, which are crucial for identifying and tracking the prevalence and dissemination dynamics of concerning viral variants. Targeted Sanger sequencing presents a fast and cost-effective means to accurately extend the coverage of whole-genome sequences. We designed a custom set of primers to amplify a 401 bp segment of the receptor-binding domain (RBD) (between positions 22698 and 23098 relative to the Wuhan-Hu-1 reference). We then designed a Sanger sequencing wet-laboratory protocol. We applied the primer set and wet-laboratory protocol to sequence 222 samples that were missing positions with key mutations K417N, E484K, and N501Y due to poor coverage after NGS sequencing. Finally, we developed SeqPatcher, a Python-based computational tool to analyse the trace files yielded by Sanger sequencing to generate consensus sequences, or take preanalysed consensus sequences in fasta format, and merge them with their corresponding whole-genome assemblies. We successfully sequenced 153 samples of 222 (69 %) using Sanger sequencing and confirmed the occurrence of key beta variant mutations (K417N, E484K, N501Y) in the S genes of 142 of 153 (93 %) samples. Additionally, one sample had the Y508F mutation and four samples the S477N. Samples with RT-PCR Ct scores ranging from 13.85 to 37.47 (mean=25.70) could be Sanger sequenced efficiently. These results show that our method and pipeline can be used to improve the quality of whole-genome assemblies produced using NGS and can be used with any pairs of the most used NGS and Sanger sequencing platforms.

SARS-CoV-2 prolonged infection during advanced HIV disease evolves extensive immune escape.

Characterizing SARS-CoV-2 evolution in specific geographies may help predict properties of the variants that come from these regions. We mapped neutralization of a SARS-CoV-2 strain that evolved over 6 months from ancestral virus in a person with advanced HIV disease in South Africa; this person was infected prior to emergence of the Beta and Delta variants. We longitudinally tracked the evolved virus and tested it against self-plasma and convalescent plasma from ancestral, Beta, and Delta infections. Early virus was similar to ancestral, but it evolved a multitude of mutations found in Omicron and other variants. It showed substantial but incomplete Pfizer BNT162b2 escape, weak neutralization by self-plasma, and despite pre-dating Delta, it also showed extensive escape of Delta infection-elicited neutralization. This example is consistent with the notion that SARS-CoV-2 evolving in individual immune-compromised hosts, including those with advanced HIV disease, may gain immune escape of vaccines and enhanced escape of Delta immunity, and this has implications for vaccine breakthrough and reinfections.

T cell responses to SARS-CoV-2 spike cross-recognize Omicron.

The SARS-CoV-2 Omicron variant (B.1.1.529) has multiple spike protein mutations1,2 that contribute to viral escape from antibody neutralization3-6 and reduce vaccine protection from infection7,8. The extent to which other components of the adaptive response such as T cells may still target Omicron and contribute to protection from severe outcomes is unknown. Here we assessed the ability of T cells to react to Omicron spike protein in participants who were vaccinated with Ad26.CoV2.S or BNT162b2, or unvaccinated convalescent COVID-19 patients (n = 70). Between 70% and 80% of the CD4+ and CD8+ T cell response to spike was maintained across study groups. Moreover, the magnitude of Omicron cross-reactive T cells was similar for Beta (B.1.351) and Delta (B.1.617.2) variants, despite Omicron harbouring considerably more mutations. In patients who were hospitalized with Omicron infections (n = 19), there were comparable T cell responses to ancestral spike, nucleocapsid and membrane proteins to those in patients hospitalized in previous waves dominated by the ancestral, Beta or Delta variants (n = 49). Thus, despite extensive mutations and reduced susceptibility to neutralizing antibodies of Omicron, the majority of T cell responses induced by vaccination or infection cross-recognize the variant. It remains to be determined whether well-preserved T cell immunity to Omicron contributes to protection from severe COVID-19 and is linked to early clinical observations from South Africa and elsewhere9-12.